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Abstract. In this paper we review the main and the most recent evidence for the 
presence of a core radius in the distribution of the dark matter around spiral galax- 
ies. Their rotation curves, coadded according to the galaxy luminosity, conform to 
an Universal profile which can be represented as the sum of an exponential thin 
disk term plus a spherical halo term with a flat density core. From dwarfs to giants, 
these halos feature a constant density region of size ro and core density po related by 
po = 4.5x 10 _2 (ro/kpc) _2//3 MQpc -3 . At the highest masses po decreases exponentially, 
with ro revealing a lack of objects with disk masses > W 11 Mq and central densities 
> 1.5 x lCP 2 (ro/kpc)~ 3 MQpc -3 , which implies a maximum mass of ~ 2 x 10 12 Mq for 
halos hosting spirals. The fine structure of dark matter halos is obtained from the kine- 
matics of a number of suitable low-luminosity disk galaxies. The inferred halo circular 
velocity increases linearly with radius out to the edge of the stellar disk, implying a 
constant dark halo density over the entire disk region. The structural properties of ha- 
los around normal spirals are similar to those around dwarf and low surface brightness 
galaxies; nevertheless they provide far more substantial evidence of the discrepancy 
between the mass distributions predicted in the Cold Dark Matter scenario and those 
actually detected around galaxies. 



1 Introduction 

Rotation curves (RC's) of disk galaxies are the best probe for dark matter (DM) 
on galactic scale. Notwithstanding the impressive amount of knowledge gathered 
in the past 20 years, only very recently we start to shed light to crucial aspects 
of the mass distribution including the actual density profile of dark halos and its 
claimed universality. 

On the cosmological side, high-resolution cosmological N-body simulations 
have shown that cold dark matter (CDM) halos achieve a specific equilibrium 
density profile [16 hereafter NFW, 6, 10, 14, 11]. This can be characterized by 
one free parameter, e.g. M200, the halo mass contained within the radius inside 
which the average over-density is 200 times the critical density of the Universe 
at the formation epoch. In their innermost region the dark matter profiles show 
some scatter around an average profile which is characterized by a power-law 
cusp p ~ r -7 , with 7 = 1 — 1.5 [16, 14, 2]. In detail, the DM density profile is: 

PNFwM = J j .,^ S j T2 (1) 

(r/r s )(l + r/r s y 
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where r s is a characteristic inner radius and p s the corresponding density. Let 
us define the halo virial radius i? v j r as the radius within which the mean density 
is Z\ v ir times the mean universal density p m at that redshift, and the associated 
virial mass M v - lr and velocity Kir = GM vlT / R vlr . By defining the concentration 
parameter as c v ; r = R v i T /r s the halo circular velocity Vcdm(^) takes the form 
[2]: 



A(x) 



where x = r/r s and A(x) = ln(l + x) — x/(l + x). As the relation between 
V v i r and i? v i r is fully specified by the background cosmology, we assume the 
currently popular /1CDM cosmological model, with fi m — 0.3, Ha = 0.7 and 
h = 0.75, in order to reduce from three to two (c v i r and r s ) the independent 
parameters characterizing the model. According to this model, Z\ v j r ~ 340 at 
z ~ 0. Let us stress that a high density fi m = 1 model, with a concentration 
parameter c v ; r > 12, is definitely unable to account for the observed galaxy 
kinematics [13]. Until recently, due to both the limited number of suitable RC's 
and to uncertainties on the exact amount of luminous matter in the innermost 
regions of spirals, it has been difficult to investigate the internal structure of 
their dark halos. However, as a result of substantial observational and theoretical 
progresses, we have recently derived the main features of their mass distribution 
for a) the Universal Rotation Curve [20] built by coadding 1000 RC's and b) a 
number of suitably selected RC's [1]. 



2 The URC and CDM Halos 



The assumed (and well supported) framework is: a) the mass in spirals is dis- 
tributed according to the Inner Baryon Dominance (IBD) regime: there is a 
characteristic transition radius Ribd — ^Rd(V op t/220 km/s) 1 ' 2 {Rd is the disk 
scale-length and V op t = V(R op t)) according which, for r < Ribd, the lumi- 
nous matter totally accounts for the mass distribution, whereas, for r > Ribd, 
DM rapidly becomes the dominant dynamical component [26, 24, 1]. Then, al- 
though the dark halo might extend down to the galaxy center, it is only for 
r > Ribd that it gives a non- negligible contribution to the circular velocity, b) 
DM is distributed in a different way with respect to any of the various baryonic 
components [20, 7], and c) HI contribution to the circular velocity at r < R op t, 
is negligible [e.g. 21]. 



2.1 Halo Density Profiles 

Reference [20] have derived from 15000 velocity measurements of 1000 RC's, 
the synthetic rotation velocities of spirals V syn (^—^, jf-), sorted by luminos- 
ity (Fig. 1, with Li the /-band luminosity and Lj/L. = 10"( M ^+ 21 - 9 )/ 5 ). Re- 
markably, individual RC's have a very small variance with respect to the corre- 
sponding synthetic curves [20, 21, 22]: spirals sweep a very narrow locus in the 
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Fig. 1. Synthetic rotation curves (filled circles with error bars) and URC (solid line) 
with its separate dark/luminous contributions (dotted line: disk; dashed line: halo). 
See [20] for details 



RC-profile/amplitude/luminosity space. On the other hand, the galaxy kine- 
matical properties significantly change with luminosity [e.g. 20], so it is natural 
to relate the mass distribution with this quantity. The whole set of synthetic 
RC's has been reproduced by means of the Universal Rotation Curve (URC) 
VuRc{ r I R op t, Li I L*) which includes: a) an exponential thin disk term [9]: 

VIurc(x) = 1.28 pV % x 2 (I K - /i#i)|i.6* (3) 
and b) a spherical halo term: 

Vl URC {x) = V 2 pt (1 -/?)(! + a 2 ) ,2 + a2) ' W 
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with x = r/Ropt, (3 = (V d ,uRc(Ropt)/V opt ) 2 , V opt = V(R opt ) and a the halo 
core radius in units of R op t- At high luminosities, the contribution from a bulge 
component has also been considered. 
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Fig. 2. a vs (3 and (3 vs V op t 



Let us stress that the halo velocity functional form (4) does not bias the mass 
model: it can equally account for maximum-disk, solid-body, no-halo, all-halo, 
CDM and core-less halo mass models. In practice, the synthetic curves V syn 
select the actual model out of the family of models Vjj RC {x) — Vfi URC (x, (3, a) + 
V 2 URC (x, (3), where a and (3 are free parameters. Adopting a ~ 1.5(i//i*) 1 / 5 
and (3 ~ 0.72+ 0.441og(L//L*) [20] or, equivalently, the corresponding a = a{(3) 
and (3 = (3(\ogV opt ) plotted in Fig. 2, the URC reproduces the synthetic curves 
V syn (r) within their r.m.s. (see Fig. 1). More in detail, at any luminosity and 
radius, \Vurc — V syn \ < 2% and the la fitting uncertainties on a and (3 are about 
20% [20]. 

To cope with this observational evidence and conveniently frame the halo 
density properties, we adopted the empirical profile proposed by Burkert [3]: 



Pb(r) 



(5) 



(r + r )(r 2 + r^) 

where po and r are free parameters which represent the central DM density and 
the scale radius. Within spherical symmetry, the mass distribution is given by: 

M b (r) = 4M {ln(l + r/r Q ) - arctan(r/r ) + 0.5 ln[l + (r/r ) 2 ]} (6) 

with Mo, the dark mass within the core, given by Mo — 1.6po r o- The halo 
contribution to the circular velocity is then: 



V b 2 (r) = GM b (r)r 



(7) 



Although the dark matter core parameters r , po and M are in principle inde- 
pendent, the observations reveal a clear correlation [3]: 
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/ \ 7/3 

M = 4.3 x !0 7 ( -rr~ ) M © ( 8 ) 

which, together with the above relationship, indicates that dark halos represent 
a 1-parameter family which is completely specified, e.g. by the core mass. 




R (kpo) 



Fig. 3. URC-halo rotation curves (filled circles with error bars) and the Burkert model 
(solid line). The bin magnitudes are also indicated 



We then compare the dark halo velocities obtained with (3) and (4), with 
the Burkert velocities Vb(r) of (5)- (7), leaving po and ro as free parameters, 
i.e. we do not impose the relationship (8). The results are shown in Fig. 3: at 
any luminosity, out to the outermost radii (~ 6Rd), Vb(r) is indistinguishable 
from Vh,uRc{ r )- More specifically, by setting Vh,uRc{ r ) = V&(r), we are able to 
reproduce the synthetic rotation curves V syn (r) at the level of their r.m.s. For 
r >> 6Rd, i.e. beyond the region described by the URC, the two velocity profiles 
progressively differ. 

The values of ro and po from the URC agree with the extrapolation at high 
masses of the scaling law p oc r^ 2 ^ 3 [3] established for objects with core radii 
ro ten times smaller (see Fig. 4). Let us notice that the core radii are very 
large: ro S> Rd so that an ever-rising halo RC cannot be excluded by the data. 

— 1/3 

Moreover, the disk-mass vs. central halo density relationship po oc M d , found 
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for dwarf galaxies [3], according to which the densest halos harbor the least 
massive disks, holds also for disk systems of stellar mass up to 10 11 M Q (see Fig. 
4). 
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Fig. 4. (up) Disk mass (in solar units) vs central halo density po (in g/cm 3 ) for normal 
spirals (filled circles) . The straight line is from [3] (bottom) central density vs core radii 
(in kpc) for normal spirals (filled circles). The straight line and the point are from the 
dwarfs sample of [3]. The curved line is: po = 5 x lCP 24 r ( ^ 2/ ' 3 exp — (ro/27) 2 g/cm 3 . 

The above relationship shows a curvature at the highest masses/lowest den- 
sities that can be related to the existence of an upper limit in the dark halo mass 
M200 Q which is evident by the sudden decline of the baryonic moss function of 
disk galaxies at M™ ax = 2 x 10 u Mq [26], that implies a maximum halo mass 
of M$§ x ~ f2 /n b M™ ax , where fi and D b ~ 0.03 [e.g. 5] are the matter and 
baryonic densities of the Universe in units of critical density. From the definition 
of M 2 oo, by means of eq. (6) and (8), we can write M 2 oo in terms of the "ob- 
servable" quantity M : Af 2 oo = ?7-&fo- For (fl$,z) — (0.3,3), ij ~ 12; notice that 
there is a mild dependence of rj on z and f2o which is irrelevant for the present 
study. From simple manipulation of previous equation- we obtain an upper limit 
for the central density, pa < 1 x 10 _20 (ro/kpc) -3 g/cm 3 , which implies a lack 

1 The virial halo mass is given by M200 = 200 x 47r/3p c i?20o^o(l + z3 )g( z ) with z the 
formation redshift, Jfeoo the virial radius, for g(z) see e.g. [2]; the critical density is 
defined as: p c = 3/(87r)G" 1 i/ 2 . 
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of objects with p > 4 x 10 25 g/cm 3 and r > 30 kpc, as is evident in Fig. 
4. Turning the argument around, the deficit of objects with Md ~ M™ ax and 
po > 4 x 10~ 25 g/cm 3 , suggests that, at this mass scale, the total-to-baryonic 
density ratio nears the cosmological value Q/Qt ~ 10. 

2.2 Testing CDM 

Out to two optical radii, the Burkert density profile reproduces, for the whole 
spiral luminosity sequence, the DM halos mass distribution. This density profile, 
though at very large radii coincides with the NFW profile, approaches a con- 
stant, finite density value at the center, in a way consistent with an isothermal 
distribution. This is in contradiction to CDM halo properties which predict [e.g. 
10] that the velocity dispersion a of the dark matter particles decreases towards 
the center to reach a — > for r — > . The dark halo inner regions, therefore, 
cannot be considered as kinematically cold structures but rather as "warm" re- 
gions with size r oc Pg ' . The halo core sizes are very large: r ~ 4 — 7Rd- 
Then, the boundary of the core region is well beyond the region where the stars 
are located and, as in [7] , even at the outermost observed radius there is not the 
slightest evidence that dark halos converge toa/)~ r~ 2 (or a steeper) regime. 

3 Individual RC's and CDM 

To derive the halo density from an individual rotation curve is certainly com- 
plicated, however, the belief according to which RC's lead to ambiguous halo 
mass modeling [e.g. 28] is incorrect. In fact this is true only for rotation curves 
of low spatial resolution, i.e. with < 3 measures per exponential disk length 
scale Rd, as for most of HI RCs. Since the parameters of the galaxy structure 
are very sensitive to the shape of the rotation curve in the region < r < Rd, 
that corresponds to the region of the RC steepest rise, then the mass model 
cannot be inferred if such a region is poorly sampled and/or radio beam-biased. 
Instead, high-quality optical RCs with tens of independent measurements in the 
critical region probe the halo mass distribution and resolve their structure. Since 
the dark component can be better traced when the disk contributes to the dy- 
namics in a modest way, it is convenient to investigate DM-dominated objects, 
like dwarf and low surface brightness (LSB) galaxies. It is well known that for 
the latter there are claims of dark matter distributions with regions of constant 
density well different from the cusped density distributions of the Cold Dark 
Matter scenario [e.g. 8, 13, 3, 4, 11, 12, 27]. However, these results are far from 
certain being 1 ) under the (unlikely) caveat that the low spatial resolution of 
the RCs does not bias the derived mass model and 2) uncertain, due to the 
limited amount of available kinematical data [see 29] . Since most of the proper- 
ties of cosmological halos are claimed universal, we concentrate on a small and 
particular sample of RCs, that, nevertheless, reveal the properties of the DM 
halos around spirals. A more useful strategy has been to investigate a number 
of high quality optical rotation curves of low luminosity late-type spirals, with 
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/-band absolute magnitudes -21.4 < Mj < -20.0 and that 100 < V opt < 170 
km s^ 1 . Objects in this luminosity/velocity range are DM dominated [e.g. 20] 
but their RC's, measured at an angular resolution of 2", have a spatial resolu- 
tion of w ~ 100(D/10 Mpc) pc and ridata ~ R op t/w independent measurements. 
For nearby galaxies: w << Rd and ridata > 25. Moreover, we select RC's of 
bulge-less systems, so that the stellar disk is the only baryonic component for 
r < R d - 

In detail, we take from [19] the rotation curves of the 'excellent' subsamplc 
of 80 galaxies, which are suitable for an accurate mass modeling. In fact, these 
RC's properly trace the gravitational potential in that: 1) data extend at least 
to the optical radius, 2) they are smooth and symmetric, 3) they have small rms, 
4) they have high spatial resolution and a homogeneous radial data coverage, i.e. 
about 30— 100 data points homogeneously distributed with radius and between 
the two arms. From this subsample we extract 9 rotation curves of low luminosity 
galaxies (5 x 1O 9 L < Lj < 2 x 1O 1O L ; 100 < V opt < 170 km s" 1 ), with their I- 
band surface luminosity being an (almost) perfect radial exponential. These two 
last criteria, not indispensable to perform the moss decomposition, are however 
required to infer the dark halo density distribution. Each RC has 7—15 velocity 
points inside R op t, each one being the average of 2 — 6 independent data. The 
RC spatial resolution is better than 1/20 R op t, the velocity r.m.s. is about 3% 
and the RC's logarithmic derivative is generally known within about 0.05. 

3.1 Halo Density Profiles 

We model the mass distribution as the sum of two components: a stellar disk and 
a spherical dark halo. By assuming centrifugal equilibrium under the action of 
the gravitational potential, the observed circular velocity can be split into these 
two components: 

V 2 (r)=V 2 (r) + V 2 (r) (9) 

By selection, the objects are bulge-less and the stellar component is distributed 
like an exponential thin disk. Light traces the mass via an assumed radially 
constant mass-to-light ratio. In the r.h.s of (9) we neglect the gas contribution 
V gas {r) since in normal spirals it is usually modest within the optical region 
[21, Fig. 4.13]: f3 gas = (V 2 as /V 2 ) Ropt ~ 0.1. Furthermore, high resolution HI 
observations show that in low luminosity spirals: V gas (r) ~ for r < Rd and 
V gas (r) ~ (2Q±5)(r-R d )/2R d for R d < r < 3R d . Thus, in the optical region: i) 
v g as ( r ) « y2 ( r ) and d[V 2 {r) - V£ as {r)]/dr > 0. This last condition implies 
that by including V gas in the r.h.s. of (9) the halo velocity profiles would result 
steeper and then the core radius in the halo density larger. Incidentally, this is not 
the case for dwarfs and LSBs: most of their kinematics is affected by the HI disk 
gravitational pull in such a way that neglecting it could bias the determination 
of the DM density. The circular velocity profile of the disk is given by (3) and 
the DM halo will have the form given by (4). Since we normalize (at R opt ) 
the velocity model (V 2 + V 2 ) 1 / 2 to the observed rotation speed V opt , (3 enters 
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explicitly in the halo velocity model and this reduces the free parameters of the 
mass model to two. 

It is important to remark that, out to R op t, the proposed Constant Den- 
sity Region (CDR) mass model of (4) is instead neutral with respect to all the 
proposed models. Indeed, by varying (3 and a, we can efficiently reproduce the 
maximum-disk, the solid-body, the no-halo, the all-halo, the CDM and the 
core-less-halo models. For instance, CDM halos with concentration parameter 
c = 5 and r s = R op t are well fit by (4) with a ~ 0.33. 

For each galaxy, we determine the values of the parameters (3 and a by means 
of a x 2 -minimization fit to the observed rotation curves: 

VLdelir; P, a) = V d 2 (r; (3) + V d 2 (r; (3, a) (10) 

A central role in discriminating among the different mass decompositions is 
played by the derivative of the velocity field dV/dr. It has been shown [e.g. 
18] that by taking into account the logarithmic gradient of the circular velocity 
field defined as: V(r) = d ^f ^"^ one can retrieve the crucial information stored 
in the shape of the rotation curve. Then, we set the x 2 -s as the sum of those 
evaluated on velocities and on logarithmic gradients: \v ~ v, ~ Vmo gy <,r "^' a ^ 

and Xv = J2"=i v(ri) ~ v ^ ei(ri:/? ' a) , with V mo del(n, (3, a) given from the above 
equations. The parameters of the mass models are finally obtained by minimizing 
the quantity Xtot = Xv + Xv- 




-0.4 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.2 0.2 0.4 

log o. 

Fig. 5. Halo parameters (a is in units of R op t) with their uncertainties 

The parameters of the best-fit models are shown in Fig. 5. They are very 
well specified: the allowed values span a small and continuous region of the (a, 
(3) space. We get a "lowest" and a "highest" halo velocity curve by subtracting 
from V(r) the maximum and the minimum disk contributions Vd(r) obtained by 
substituting in (3) the parameter (3 with f3b es t + 5(3 and [3b es t — §0, respectively. 
The derived mass models are shown in Fig. 6, alongside with the separate disk 
and halo contributions. It is then obvious that the halo curve is steadily increas- 
ing, almost linearly, out to the last data point. The disk-contribution (3 and the 
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halo core radius a span a range from 0.1 to 0.5 and from 0.8 to 2.5, respectively. 
In each object the uniqueness of the resulting halo velocity model can be realized 
by the fact that the maximum-disk and minimum-disk models almost coincide. 
Remarkably, we find that the size of the halo density core is always greater than 
the disk characteristic scale-length Rd and it can extend beyond the disk edge 
(and the region investigated). 




Fig. 6. CDR model fits {thick solid line) to the RCs {points with errorbars). Thin solid 
lines represent the disk and halo contributions. The maximum disk and the minimum 
disk solutions are also plotted {dashed lines) 



3.2 Testing CDM 

In Fig. 7 we show the halo velocity profiles for the nine galaxies. The halo circular 
velocities are normalized to their values at R op t and expressed as a function of 
the normalized radius r/R opt . These normalizations allow a meaningful compar- 
ison between halos of different masses. It is then evident that the halo circular 
velocity, in every galaxy, rises almost linearly with radius, at least out to the 
disk edge: V h (r) cx r for 0.05i? opt < r < R opt . 

The halo density profile has a well defined (core) radius within which the 
density is approximately constant. This is inconsistent with the singular halo 
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density distribution emerging in the Cold Dark Matter (CDM) halo formation 
scenario. More precisely, since the CDM halos are, at small radii, likely more 
cuspy than the NFW profile: pcdm oc r -1 - 5 [e.g. 14], the steepest CDM halo 
velocity profile Vh(r) oc r 1 / 4 results too shallow with respect to observations. 
Although the mass models of (4) converge to a distribution with an inner core 
rather than with a central spike, it is worth, given the importance of such result, 
also checking in a direct way the (in)compatibility of the CDM models with 
galaxy kinematics. We assume the NFW two-parameters functional form for the 
halo density [15, 16, 17], given by (1). Though N-body simulations and semi- 
analytic investigations indicate that the two parameters c v i r and r s correlate, 
they are left independent to increase the chance of a good fit. For the object 
under study a generous halo mass M V i r upper limit is M up = 2 x 1O 12 M . 




0.2 0.4 0.6 o.a 



r/R opl 

Fig. 7. The halo velocity profiles of the sample spirals. Vh(r) rises almost linearly with 
radius: the DM halo density remains approximately constant 

The fits to the data are shown in Fig. 8 and compared with the NFW mod- 
els: for seven out of nine objects the latter are unacceptably worse than the 
CDR solutions, moreover in all objects, the CDM virial mass is too high high: 
M V i r ~ 2 x 10 12 M Q and the resulting disk mass-to-light ratio too low. The in- 
adequacy of the CDM model for our sample galaxies is even more evident if one 
performs the fit after removing the constraint on virial mass. In fact, good fits are 
obtained only for very low values of the concentration parameter (c v i r ~ 2) and 
for ridiculously large virial velocities and masses {V v -„ ~ 600—800 km s _1 ; M vir ~ 
10 13 — 10 14 M Q ). These results can be explained as effect of the attempt, by the 
minimization routine, to fit the NFW velocity profile (V(r) oc r 5 ) to data 
intrinsically linear in r. 

4 Conclusions: an Intriguing Evidence 

The dark halos around spirals emerge as an one-parameter family; it is relevant 
that the order parameter (either the central density or the core radius) corre- 
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lates with the luminous mass. However, we do not know how it is related to the 
global structural properties of the dark halo, like the virial radius or the virial 
mass. The halo RC, out to 6Rd, is completely determined by parameters, i.e. 
the central core density and the core radius, which are not defined in present 
gravitational instability/ hierarchical clustering scenario. In fact the location of 
spiral galaxies in the parameter space of virial mass, halo central density and 
baryonic mass, determined by different processes on different scales , degenerates 
with no doubt into a single curve (see Fig. 4), we recall that: Po = -^20o/ r o 
and M d — G^ 1 (3V^ pt R op t, of difficult interpretation within the standard theory 
of galaxy formation. Crucial insight has come from disk-halo density decompo- 




Fig. 8. NFW best-fits solid lines of the rotation curves (filled circles) compared with 
the CDR fits (dashed lines). The x values are also indicated 



sitions of a number of disk galaxies. These galaxies have a relevant amount of 
dark matter: the contribution of the luminous matter to the dynamics is small 
and it can be easily taken into account. Moreover, the high spatial resolution 
of the available rotation curves allows us to derive the separate dark and lumi- 
nous density profiles. We find that dark matter halos have a constant central 
density region whose size exceeds the stellar disk length-scale Rd- As result, the 
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halo profiles disagree with the cuspy density distributions typical of CDM halos 
which, therefore, fail to account for the actual DM velocity data. 

Pointing out that a review on the various efforts aimed to cope with the core 
radii evidence will be published elsewhere, we conclude by stressing that, for any 
theory of galaxy formation, time is come to seriously consider that stellar disks 
(and perhaps also stellar spheroids) lay down in dark halos of constant density. 
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